DEPINNING TRANSITIONS IN DISCRETE REACTION-DIFFUSION 

EQUATIONS * 

A. CARPIOt AND L. L. BONILLAtS 

Abstract. We consider spatially discrete bistable reaction-diffusion equations that admit wave 
front solutions. Depending on the parameters involved, such wave fronts appear to be pinned or 
to glide at a certain speed. We study the transition of traveling waves to steady solutions near 
threshold and give conditions for front pinning (propagation failure). The critical parameter values 
are characterized at the depinning transition and an approximation for the front speed just beyond 
threshold is given. 
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1. Introduction. Spatially discrete systems describe physical reality in many 
different fields: atoms adsorbed on a periodic substrate |l3|, motion of dislocations 
in crystals propagation of cracks in a brittle material [35|, microscopic theories 
of friction between solid bodies |Q , propagation of nerve impulses along myelinated 
fibers |2^, pulse propagation through cardiac cells ||2^, calcium release waves 
in living cells 0, sliding of charge density waves superconductor Josephson 

array junctions ^9|, or weakly coupled semiconductor superlattices No one 

really knows why, but spatially discrete systems of equations often times have smooth 
solutions of the form Un(t) — u{n — ct), which are monotone functions approaching 
two different constants as (n — ct) —> ±oo. Existence of such wave front solutions 
has been proved for particular discrete systems having dissipative dynamics po|] . In 
the case of discrete systems with conservative dynamics, a wave front solution was 
explicitly constructed by Flach et al. A general proof of wave front existence for 
discrete conservative systems with bistable sources is lacking. 

A distinctive feature of spatially discrete reaction-diffusion systems (not shared 
by continuous ones) is the phenomenon of wave front pinning: for values of a control 
parameter in a certain interval, wave fronts joining two different constant states fail 
to propagate ||2^. When the control parameter surpasses a threshold, the wave front 
depins and starts moving ||2^, |l^, 32 , |j . The existence of such thresholds is thought 



to be an intrinsically discrete fact, which is lost in continuum aproximations. The 
characterization of propagation failure and front depinning in discrete systems is thus 
an important problem, which is not yet well understood despite the numerous inroads 
made in the literature H ||, ||, ||, H, |2|, |3|, H, |3|. 

In this paper, we study front depinning for infinite one-dimensional nonlinear 
spatially discrete reaction-diffusion (RD) systems. When confronted with a spatially 
discrete RD system, a possible strategy is to approximate it by a continuous RD 
system. For generic nonlinearities, the width of the pinning interval is exponentially 
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small as the continuum limit is approached. Pinning in the continuum limit has been 
analyzed by many authors using exponential asymptotics, also known as asymptotic 
beyond all orders. As far as we can tell, usage of these techniques for discrete equations 
goes back to two classic papers by V.L. Indenbom |^ (for the FK potential) and by 
J.W. Cahn [Q (for the double- well potential). In both cases, an exponential formula 
for the critical field was derived by means of the Poisson sum rule. In the context of 
dislocation motion, exponential formulas for the depinning shear stress of the Peierls- 
Nabarro (PN) model were found earlier by Peierls and Nabarro . Descriptions 
of wave front pinning near the continuum limit can also be found in more recent work 

©11,13; 

Analyzing the continuum limit of a discrete system by means of exponential 
asymptotics is a costly strategy to describe pinning for two reasons. It is not nu- 
merically accurate as we move away from the continuum limit, and it ceases to be 
useful if convective terms ||ri| or disorder ||l2| alter the structure of the discrete sys- 
tem (quite common in applications). Thus other authors have tried to describe the 
opposite strongly discrete limit. For discrete RD equations, Erneux and Nicolis 
studied a finite discrete RD equation with a cubic nonlinearity, Dirichlet boundary 
condition at one end and a Neumann boundary condition at the other end. They 
considered a particular limit in which two of the three zeroes of the cubic nonlinearity 
coalesced as diffusivity went to zero. Erneux and Nicolis's calculation is essentially a 
particular case of our active point approximation that involves only one active point 
and makes an additional assumption on the nonlinearity (not needed in our calcula- 
tions). They found that the wave front velocity scales as the square root of (d — dc) 
{d is the diffusivity and dc its critical value at which wave fronts are pinned). Es- 
sentially the same results can be found in the Appendix of Ref. Kladko et al 
Ipsf introduced an approximation called the single active site theory. In this approx- 
imation, the wave front is described by two linear tails (solution of the RD equation 
linearized about each of the two constants joined by the front) patched at one point. 
This approximation was used to estimate the critical field for wave front depinning. 

By a combination of numerical and asymptotic calculations, we arrive at the 
following description The nature of the depinning transition depends on 

the nonlinearity of the model, and is best understood as propagation failure of the 
traveling front. Usually, but not always, the wave front profiles become less smooth 
as a parameter F (external field) decreases. They become discontinuous at a critical 
value Fc- Below Fj., the front is pinned at discrete positions corresponding to a stable 
steady state. As a consequence of the maximum principle for spatially discretized 
parabolic equations, stationary and moving wave fronts cannot simultaneously exist 
for the same value of ||^ . This is not the case for chains with conservative dynamics, 
that are spatially discretized hyperbolic equations without a maximum principle. For 
chains with conservative Hamiltonian dynamics, an inverse method due to Flach, 
Zolotaryuk and Kladko | p6| explicitly shows that stationary and moving wave fronts 
may coexist for the same value of the parameters. 

We consider chains of diffusively coupled overdamped oscillators in a potential F , 
subject to a constant external force F: 

du 

(1.1) — ^ = Un+l - 2Un + Un-1 + F - Ag{Un). 

dt 



Here g{u) = V'{u) is at least and it presents a 'cubic' nonlinearity (see Fig. LI), 
such that Ag{u) - F has three zeros, Ui{F/A) < U2{F/A) < Uz{F/A) in a certain 
force interval {g'{U^{F/A)) > for i = 1,3, g'{U2{F/A)) < 0). Provided g{u) is odd 
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with respect to C/2(0), there is a symmetric interval \F\ < Fc where the discrete wave 
fronts joining the stable zeros Ui{F/A) and UsiF/A) are pinned j2|, ^. For \F\ > Fc, 
there are smooth travehng wave fronts, u„(t) = u{n — ct), with u(— oo) = Ui and 
u{oo) = Uz, as proved in The velocity c{A,F) depends on A and F and 

it satisfies cF < and c ^ as |F| ^ Fc M. Examples are the overdamped 



Frenkel-Kontorova (FK) model {g = sinu; see Fig. LI) ||17[] and the quartic double 
well potential (V — (u^ — 1)^/4). Less symmetric nonlinearities yield a non-symmetric 
pinning interval and our analysis applies to them with trivial modifications. Note 
that coexistence of fronts traveling in opposite directions can occur in the case of 
conservative systems, but not for Eq. (1.1) due to the maximum principle (which is 
the basis of comparison techniques) 



For the overdamped FK model given by Eq. (LI) with g — sinu, Fig. L2 shows 
wave front profiles near the critical field. Individual points undergo abrupt jumps at 
particular times, which gives the misleading impression that the motion of the dis- 
crete fronts proceeds by successive jumps. Actually, the points remain very close to 
their stationary values at = Fc, say Un{A,Fc), during a very long time interval of 
order \F — Fc\~2. Then, at a specific time, all the points u„(t) jump to a vicinity of 
Un+i{A, Fc). The method of matched asymptotic expansions can be used to describe 
this two-stage motion of the points Un{t). Then the wave front profile can be recon- 
structed by using the definition M„(i) = u{n — ct). The slow stage of front motion is 
described by the normal form of a saddle-node bifurcation and it yields an approxi- 
mation to the wave front velocity, that scales with the field as \F — F^^ . This scaling 
has been mentioned by other authors: it was found numerically in 0| and by means of 
exponential asymptotics in the limit A small in p7[ |. It is also conjectured in on 
the correct basis that the depinning transition consists of a saddle-node bifurcation 
(a similar claim was stated in pO[ | for continuous reaction diffusion equations with 
localized sources). However, the derivation of the local saddle-node normal form and 
the correct description of the global saddle-node bifurcation involving matching with 
a fast stage during which the front jumps abruptly one lattice period were apparently 
omitted by the authors of Ref. who used energy arguments. Our picture of the 
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Fig. 1.2. Wave front profiles for the overdamped FK model when A = 1 near Fc, 



wave front depinning transition has essentially been corroborated by King and Chap- 
man in the continuum limit (as an appropriate dimensionless lattice length goes to 
zero) who used asymptotics beyond all orders An independent confirmation fol- 
lows from Fath's calculations for a spatially discrete reaction-diffusion equation with 
a piecewise linear source term ||l5[ (except that the velocity should scale differently 
with \F — Fc\ in this case). 

For exceptional nonlinearities, the wave front does not lose continuity as the field 
decreases. In this case, there is a continuous transition between wave fronts moving to 
the left for F > and moving to the right for F < 0: as for continuous systems, front 
pinning occurs only at a single field value F — |l6|, ^ . Wave front velocity 

scales then linearly with the field. We discuss the characterization of the critical 
field (including analytical formulas in the strongly discrete limit), describe depinning 
anomalies (discrete systems having zero critical field ^), and give a 

precise characterization of stationary and moving fronts near depinning (including 
front velocity) by singular perturbation methods. Our approximations show excellent 
agreement with numerical simulations. 

The rest of the paper is organized as follows. In section ||, we characterize wave 
front depinning. We also explain that pinning of wave fronts normally occurs at force 
values belonging to an interval with nonzero length. However, there are nonlinearities 
for which pinning occurs only at F = Fc = 0. In Section we present a theory of 
wave front depinning for the strongly discrete case {A large). This theory enables us 
to predict the critical field and the speed and shape of the wave fronts near thresh- 
old. The main ideas of our theory are very simple. First of all, a wave front profile 
Unit) = u{n — ct) can be reconstructed if we follow the mo tion of one point during a 
sufficiently long time interval. Secondly, the analysis of Eq. (1.1) is complicated by the 
presence of the discrete diffusion term, Wn+i — 2w„-|-w„_i. Previous authors have tried 
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to approximate this term by its continumn limit (corresponding to A ^ 0), which 
leads to using exponential asymptotics (see also |^ on using exponential asymp- 
totics for the Hamiltonian version of our model). However, we are only interested 
in constructing solutions of Eq. ( |l . l| ) joining constant values. For sufficiently large 
A (say ^ = 0.1 for the FK model), Ui is approximately either Ui{F/A) or {/^{F/A) 
except for a finite number of points (the active points). Then we can approximate 
the infinite system (1.1) by a closed system of ordinary differential equations (only 
one equation for A > 10 in the FK model). The depinning transition is a global 
bifurcation of this system, as explained in Section ^. Some auxiliary technical results 
are collected in Appendix A. 

2. Fron t p inning as propagation failure. To describe monotone stationary 
solutions of (|li1 ) joining Ui{F/A) and UsiF/A) for |F| < Fc, it is better to start by 
considering traveling wave fronts for j.F'l > F^. It has been proved (and corroborated 
by numerical calculations) that traveling wave fronts and stationary profiles cannot 
coexist at the same value of F Furthermore, numerical computations of wave 
fronts near the critical fields Fc for the FK and other usual potentials show staircase- 
like wave front profiles, which sharpen as F approaches Fc- At F = Fc, a series of 
gaps open up and one is left with a discontinuous stationary profile s{x) solving 



(2.1) 



s{x + 1) - 2s{x) + s{x - 1) = Ag{s{x)) - Fc, x e H, 
s(-c») = t/i(FcM), s(c») = C/3(FcM). 



The profile s{x) is increasing and piecewise constant. The sequence of constant values 
attained by s{x) defines a steady solution u„ of (1.1) with F = Fc- A stationary 
solution can thus be understood as a wave front that fails to propagate and is pinned 
at discrete values. Fig. L2 illustrates the pinning transition for the FK model with 
A = 1. As F decreases from 0.02 to 0.0127, a series of steps are formed. Fig. L2 (c) 



depicts the paths described by three consecutive points. All profiles look identical and 
are obtained by shifting anyone of them some multiple of a certain constant length. 
This implies that the length of all steps in the profile is the same and all the points 
Un(t) in (LI) proceed to climb the next step in the staircase at the same time. This 
behavior indicates that the wave front is a traveling wave, M„(i) = u{n — ct). Proofs 
of this fact for some sources can be found in Ref. |Q . 

2.1. Limiting front profile at the critical field. Let us start by showing that 
the limit of the traveling waves as F — > Fc is singular if Fc > 0. This fact can be 
guessed from the differential-difference equations satisfied by the wave profiles. The 
traveling waves for |F| > Fc have the form u„(i) = u{n — ct), where the profile u{z) 
solves H 



(2.2) 



cuz = u{z + 1) - 2u{z) + u{z - 1) - Ag{u{z)) + F z eH, 
u{^^)^Ui{F/A), u{^) = U:,{F/A). 



The solution u is as smooth as allowed by g{u) {u is (7'^+^ if g{u) is C*^, with A: > 1). 
Mutiplying then Eq. (2.2) by Uz and integrating it, we get 



(2.3) 



tdz = F 



A first obvious conclusion is that the sign of c is opposite to the sign of F. Let 
Fc be positive. As F -> Fc, c -> and F [UsiF/A) - Ui{F/A)] Fc [UsiFc/A) - 
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Ui{Fc/A)] 7^ 0. Therefore the integrals J u^dz ^ oo as — > F^- Thus the limiting 
profile must be discontinuous if -Fc > 0. 



If Fc = 0, the relation (2.3) can be used to show that Eq. (2.1) has a smooth 



solution. In fact, provided c ~ —KF {K > 0) as F ^ 0, we can use (US) to 



bound uniformly the derivatives of the solutions u in (2.2) for F 0. Then we 
obtain a smooth solution of Eq. ( |2.lD in the limit as F ^ 0. We will come back 
to this question later on in Subsection 2.3. Note that the stationary equation s{x + 
1) — 2s{x) + s{x — 1) = Ag{s) — F has no continuous solutions joining Ui{F/A) to 
UslF/A) unless F = 0. To see this we multiply the equation by Sx (in the sense 
of distributions if necessary) and integrate to get F = 0. 

2.2. Characterization of the critical field. Some results are available in the 
continuum limit A — > 0. For g — sinw, it is well known that Fc vanishes exponentially 
fast as A goes to zero. An exponential formula for F^ was first found by Indenbom 
p2t using the Poisson sum rule (following the calculations of the PN energy barrier for 
the PN model by Peierls and Nabarro |3^), and numerically checked by Hobart 
PH in the context of the Peierls stress and energy for dislocations. For the discrete 
bistable RD equation, Cahn j^] derived an exponential dependence of Fc by a similar 
technique. Related ideas can be found in Kladko et al |E8|. These arguments can be 
used for other potentials and suggest that F^ ^ C e~''/^ as A ^ 0+ (with positive 
C and 77 independent of A) holds for a large class of nonlinearities. Using exponential 
asymptotics, King and Chapman ||2^ have obtained precise formulas for the critical 
field and the wave front velocity of a discrete RD equation. Particularized to the 
FK potential, their formulas for the critical field and for the wave front velocity after 
depinning are F^ - Ae-'''/[2«'"h-i(VA/2)]^ ^ ^ 355 ^^^^ ^ _ D {F"^ ~ F"^) / A, 

respectively. This later result agrees with the scaling law c ^ \F — Fc\^ ^ found in a 
large class of discrete RD equations |l^, ^ ^ and in continuous equations with 
localized sources |Q. However, exponential asymptotics does not work for A 
large. We shall therefore follow a different approach. We shall begin by considering 
stationary increasing discrete front profiles and study under which conditions they 
start moving. Since stationary fronts are pinned wave fronts, we can call the transition 
from stationary to moving fronts the depinning transition. 

Two facts distinguish the depinning transition: (i) the smallest eigenvalue of 



(1.1) linearized about a stable stationary profile becomes zero (see below), and (ii) 



stationary and moving wave fronts cannot coexist for the same values of the field. 



First of all, the following comparison principle |2^ for (IT) can be used to show that 
stationary and traveling wave fronts cannot coexist for the same value of F : 

Comparison principle: Assume that we have two configurations Wn{t) and 
ln(t). If initially Wn(0) > ^n(O), Vrt and at any later time t>0 

(2.4) > Wn+i - 2Wn + Wn-l - Ag{Wn) + F, 

dt 

(2.5) ^ < In+l - 2ln + In-l - Ag{ln) + F, 

dt 



for all n G then necessarily w„(t) > ^„(i) for all n and t. Wn satisfying (2.4) is 



said to be a .supers olution. Z„ satisfying (2.5) is said to be a subsolution. 



Front pinning can be proved using stationary sub and supersolutions, which can be 
constructed provided the stationary solution is linearly stable. The smallest eigenvalue 



of the linearization of (1.1) about a stationary profile Un{A,F), u„(t) = Un{A, F) 



Depinning in discrete equations 



7 




Fig. 2.1. (a) Critical field, Fc{A), for A g ^^,10^. We have compared the result of approx- 
imately solving \i{A,F) = for F as a function of A (see Appendix A) to the asymptotic result 
Fc ~ 356.1 e-'^^/[2sinh-i(v^/2)] of Ref. (b) Relative error of the exponential asymptotics 

approximation. 



VnC , is given by 

E[ivn+i-v„)^ + Ag'MA,F))vl 



(2.6) Xi{A,F) =MiTi^ 



over a set of functions u„ which decay exponentiaUy as n — > ±oo. We show in the 
Appendix that the minimum is attained at a positive eigenfunction. 

The critical field can be uniquely characterized by Xi{A, Fc) ~ and Xi{A, F) > 
for \F\ < Fc- The details are given in Appendix A. Notice that Xi{A, F) > implies 



that Eq. (2.1) does not have smooth solutions s{x). For otherwise, u„ — s'{n) is an 
eigenfunction corresponding to Ai = as it happens in the continuum limit. The 
previous characterization is the base of a procedure to calculate Fc{A). In Section 
||, we shall show that wave fronts near the depinning transition are described by a 
reduced system of equations for a finite number of points it„(i) which "jump" from 
about a discrete value corresponding to the stationary solution, m„(A, Fc), to the next 
one, Un+i{A,Fc), during front motion. The smallest eigenvalue for the linearization 
of the reduced system of equations about a stationary solution approximates well Ai. 



The critical field obtained by this procedure has been depicted in Fig. for the 
FK potential and compared to King and Chapman's asymptotic result (obtained by 
keeping two terms in their formulas). Notice that the asymptotic result loses accuracy 
as A increases. 



Equation (2.6) shows that the critical field is positive for large A and typical 
nonlinearities. In fact, consider the FK potential. For F ~ there are two one- 
parameter families of stationary solutions which are symmetric with respect to U2 
(see Fig. p^ ), one taking on the value U2 (unstable dislocation), and the other 
one having w„ ^ C/2 (stable dislocation) ||2^, ^. The centers of two stable (or two 
unstable) dislocations differ in an integer number of lattice periods. Except for a 
possible rigid shift, the stable dislocation, Un{A, 0), is a dynamic ally stable stationary 



solution towards which step-like initial conditions evolve. Figs. |2.3| (a) and (b) show 
two initial conditions that evolve (exponentially fast) towards the stable dislocation. 
Half the initial points, u„(0), have been selected to be below U2, and the other half 



are above this value. In Fig. 2.3(a), u„(0) — u„(A, 0) = e„, where e„ are real random 
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Fig. 2.2. Stable and unstable dislocations for the FK model when F = and A = 0.1. 




-10 



(a) F=0 



OOOOOOOOO 



-O- 

ooooooooo 
^ ^ ^ ^ * , 



4 



2 

-10 



(c) F=6.1 




n 




n 



10 



10 



6 
4 



(b) F=0 



* * * * 



,000000000. 



OOOOOOOOO 

-2 



-10 

8~ 
6 
4 




n 



10 



(d) F=6.1 



^0000000000 



2 WW W i 

c 

oooooooo 



-10 




n 



10 



Fig. 2.3. Initial condition Mn(0) (asterisks) and its large time limit, the stable dislocation 
(circles), for the FK model with A = 10. F = for (a) and (b) and F = 6.1 < Fc for (c) and (d). 
The initial points are selected as indicated in the text. 



numbers with |e„| < 0.5. In Fig. 2.3(b), u„(0) 
and Sn takes randomly on the values 1 or -1. 



- C/1,3 = SnB, 0<B = U2-Ui- 0.2, 
By using comparison methods, it is 



possible to prove that a small disturbance of the stable dislocation evolves towards it. 
The same results hold for the stable stationary solution w„(A, F) for < |i^| < Fc. As 
\F\ increases, a disturbance of the stable stationary solution typically evolves towards 
the same stationary solution displaced an i nte ger number of lattice periods unless the 
disturbance is sufficiently small. See Fig. 2_^(d) for an example of this phenomenon 
for F sl ight ly smaller than Fc. Carefully selecting the initial condition avoids this, as 
in Fig. |]|(c). 

For large A, the stable dislocation has g'{un) > for all n, and (|2.6|) gives 
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Fig. 2.4. Stationary solutions for the FK model when = 1, 2, 10, 100. 



Ai(^, 0) > 0. Since Ai(A, fe) = 0, this implies that the critical field is nonzero. 
(Different proofs are given in |2^, H). As A > decreases, several u„ may enter the 
region of negative slope g'(u): the number of points with g'{un) < increases as A 



decreases, see Figures p.2| and 2.4 . It should then be possible to have Ai(A, 0) — 0, i.e. 



Fc — 0, for a discrete system! Examples of this pinning anomaly will be given next. 

2.3. Pinning Failure. Despite widespread belief, it is not true that the critical 
field is positive for all discrete systems. This point was already raised by Hobart , 
who proposed the following numerical criterion to check whether for a given source g 



the critical field for (1.1) is zero. 

Let us assume for the sake of simplicity that g is odd about 0. Then ?72(0) = 
and Ui{Q) — —1/3(0). For any x G ([/i(0), C/3(0)), we can compute numerically a 
unique value y{x) such that the sequence w„ defined by uo = x, ui = y{x) and 
Un = 2u„_i — Un-2 + g(un-i),n > 1 tcuds to U^iO) a,s n ^ oo. Hobart conjectured 
that Fc = for a given nonlinearity g, provided that the function y{x) satisfies: 

(2.7) y^^{x) = ~y{-x), y{x) - y{-x) ^ 2x + g{x), 

for X € {U I (0), 1/3(0)). It is fairly easy to construct examples of nonlinearities g{x) 
for which ( [2.7] ) holds. It suffices to choose some smooth odd increasing function 
u{x) such that u{x) — ^ ±a as x — > ±00 for some a > 0. We define g{u{x)) = 
u{x + 1) — 2u{x) + u{x — 1) so that g{z) — u{u^^{z) + 1) — 2z + u{u~^{z) — 1) and 
y{z) — u{u^^{z) + 1). Choosing u{x) — tanh(a;) ||; we get an explicit formula 
for g: g{z) = —2^z{l — z^) / {1 — jz^) with 7 = tanh (1). Notice that one or two 
points of the stationary solu tion s, Un — tanh(n + p) {p is any constant), enter the 
region where g' < 0; see Fig. |2.5| (a). 

By following this procedure, we find examples of bistable source terms for which 



(1.1) has a uniparametric family of continuous stationary solutions, u„ = u{n + p), 
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< p < a, satisfying — 2u„ + — g{un) and w_oo — —a, u^o = a. In this 



case, (1.1) does not have stationary solutions joining Ui[F/A) and U^{F/A) unless 
F — ^1 . The existence of continuous steady solutions for = implies that there 
is a continuous transition from wave fronts traveling to the left (c < 0) for > to 
wave fronts traveling to the right (c > 0) for F < 0. Only at F = wave fronts are 
stationary (pinned). This pinning anomaly is stated more precisely as: 

Theorem. Let g ^ he as in the Introduction with g{0) = 0, and let C{F) be 
the operator 

(2.8) C{F) Vn = Ag'{Un)Vn + 2Vn - Vn+1 - Wn-1, 



corresponding to the evolution equation (1.1) linearized about the stationary solution 



Un = u„(A, F) at field F. Let us assume that for F = 0, there exists a differentiable 
increasing stationary solution u{x), such that u{x) iU^iO) as x ±oo. Then, 
1. Zero is the s mal lest eigenvalue of the operator Cq — C{0), corresponding to the 



evolution equation (1.1) linearized about the stationary solution u„(A, 0) — u{n). 

2. F,{A) = for ilj). 

3. Traveling wave fronts exist for all F ^ 0. Furthermore, their speed increases 
linearly with the force for small F. We have 

(2.9) e^-F^Mt^, 
as F -^0. 

Mo reover, Statement 3 implies the existence of steady differentiable solutions u{x) 



of ( TJ ) such that u{x) ±1/3(0) as x ^ ±cxd for F — 0. 

It is not our goal here to give a rigorous proof of this result, but to sketch the main 
ideas. First of all, note that the derivative Vn — Ux{n) > is a positive eigenfunc- 
tion of the elliptic operator Cq corresponding to the eigenvalue A = and decaying 
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exponentially at infinity. Statement 1 immediately follows. This fact can be used to 



construct propagating sub and supersolutions for (1.1) which forbid pinning for any 
F ^ 0. Thus, Fc — 0, which is Statement 2. For F ^ e > sufficiently small, the 
propagating subsolutions are ln(t) = /(n + ecgt), with cq > and l{x) = u{x) + eux{x). 
For F = —€, the propagating supersolutions are Wn{t) — w{n — ecot) with cq > and 
'w{x) = u{x) — eux{x). In both cases, we have to choose co < l/niax{ux)- A subsolu- 
tion traveling to the left "pushes" the fronts to the left. Similarly, the supersolutions 
traveling to the right "push" the fronts to the right. 

Let us now obtain Statement 3. If \F\ > 0, we have traveling wave front solu- 



tions Un{t) —U{n — ct) of (IT), whose profile Z//(z) satisfies the differential-diflference 
equation: 

(2.10) -c^(z) ^U{z + l) - 2U{z)+U{z - I) - g{U{z)) + F, 

dz 

and ^^(zboo) = ±Uz{Q), §]. Let F = Foe with < e < 1. The travehng wave solution 
can be written aslA{n — ct) ~ u{n — ct) + ew{n — ct) + o(e), where u(x) is the smooth 
stationary profile. Let z — n ~ ct and c = cqe + o(e). Then w obeys 

du 

w{z + 1) - 2w{z) + w{z - 1) - Ag' {u{z))w{z) = -cq — (z) - Fq, 

dz 

w{-oo) = w{oo) = ^ = ^ . 

g'(Ui) g'{U3) 

By the Fredholm alternative, this linear nonhomogeneous equation has a solution if 
the left hand side —codu/dz — Fq is orthogonal to the eigenfunction du/dz, which 



yields ( p.9D 



In Section |3|, w e show that the wave front speed c scales as \F—Fc\^ if Fc > 0. Our 
linear scaling (|2.9|' ) of the velocity in Statement 3 therefore implies that Fc = 0. The 
linear scaling ( ^.9[ ) with Fc = implies the existence of smooth stationary solutions 
at F = as discussed in the first Subsection. 



Remark 1. We conjecture the three statements in the above theorem are equivalent. 
To prove this, it would be enough to show that Fc = implies the linear scaling of 
the speed of the waves (Statement 3). Then, existence of diffcrentiable stationary 
solutions follows. This implies statement 1 (Ai(A,0) = 0), which implies Statement 
2 (Fc = 0), as we showed above. 

Remark 2. When stationary wave front solutions have smooth profiles, pinning failure 
occurs for discrete RD equations and for discrete equations with conservative dynam- 
ics. In the latter case, translation invariant smooth profiles have the same energy and 
therefore the PN energy barrier (defined as the smallest energy barrier to overcome for 
a kink or wave front to move [Q) vanishes. Pinning of a wave front usually resu lts i f 
the energy difference between the stable and the unstable front solutions (see Fig. |2!^ ) 
is not zero. This energy difference provides an estimation of the PN energy barrier. 
Discusions of the PN potential and the PN barrier can be found in Section 2.3 of 
Ref.|^ and in Section III.B of Ref. [|l^. The mathematical meaning and usefulness of 
the PN barrier for an infinite system with conservative dynamics are worth studying. 



Remark 3. Speight and Ward |^ and Speight |3|] have developed a technique 
to discretize some continuum conservative models in such a way that kink-like initial 
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A=10,F=4 

5- 



^ CT i : i m i mm# 
& 



6 (b) 

A-10,F=6.102 

5 



Fig. 3.1. Stationary solutions for the FK model with A = 10.' (a) No points are found in the 
region g' < for sufficiently small F; (b) one point enters the region g' < for sufficiently large 
F < Fc. 



profiles may propagate without getting trapped. Their idea is to seek a discrete 
version of the potential energy which admits minimals satisfying a first order difference 
equation called the Bogomol'nyi equation so that there is no PN barrier. On the other 
hand, the difference operators in Speight et al discretized equations of motion have a 
structure different from discrete diffusion and are hard to justify physically. 

Remark 4- In discrete RD equations, moving and pinned fronts cannot coexist for the 
same value of the applied field. For chains with conservative Hamiltonian dynamics, 
the situation is less clear. In fact, it is possible to have two stationary front solutions 
with a positive energy difference between them (which would imply a nonzero critical 
field and therefore wave front pinning according to general belief), and yet a moving 
wave front may coexist with the stationary fronts for the same parameter values. 
An explicit example of this situation has been constructed by Flach, Zolotaryuk and 
Kladko using an inverse method |l^ . 

3. Asymptotic theory of wave front depinning. In this section we introduce 
a systematic procedure to derive analytic expressions for the critical field > as a 
function of A, and for the front profiles and their velocity as functions oi F — Fc and 
A. Our methods work best in the strongly discrete case, for large A. Our ideas are 
quite general and may be applied successfully to more complex discrete models [ pj] . 
We shall assume that g £ throughout this Section. 

3.1. Theory with a single active point. We choose A large enough for the 



stable dislocation in Fig. not to enter the region where g' < 0, see Fig. 2.4. When 
F > 0, this solution is no longer symmetric with respect to U2. If F is not too large, 
all Un{A,F) avoid the region of negative slope g'{u) < 0. For larger F and generic 
potentials (FK, double- well,. . . ), we have observed numerically that g' < fo r a s ingle 



point, labelled uo{A,F). This property persists until Fc is reached, see Fig. B.l, 

Firstly, consider the symmetric stationary profile with u„ ^ U2 for F = 0. The 
front profile consists of two tails with points very close to Ui and U3, plus two sym- 
metric points uq, Ml in the gap region between Ui and U3. As F > increases, this 
profile changes slightly: the two tails are still very close to Ui{F/A) and U3{F/A). 
As for the two middle points, ui gets closer and closer to U3 whereas uq moves away 
from Ui. This structure is preserved by the traveling fronts above the critical field: 
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there is only one active point most of the time, which we can adopt as our uq. Then 
the wave front profile (2^) can be calculated as u{—ct) = UQ{t). In Eq. (pTl]), we can 
approximate u_i ^ J/i, ui ~ L/3, thereby obtaining 



(3.1) 



duo 
IT 



Us 



2uo-Ag{uQ)+F. 



This equation has three stationary solutions for F < F^, two stable and one unstable, 
and only one stable stationary solution for F > Fc- Let us consider F < Fc- Only 



two out of the three solutions of Eq. (3.1) approximate stationary fronts for the exact 



system: those having smaller values of uq. The one having smallest uq approximates 
the stable stationary front, the other one approximates the unstable stationary front. 
Recall that the unstable front had a value uq = [Ui{0) + J73(0)]/2 at the middle of 
the gap for F = 0. As > increases, mq decreases towards Ui{F/A). Thus one 
active point will also approximate the profile of the unstable stationary front. The 
stationary solution of Eq. (3.1) having the largest value of uq (slightly below U3{F/A)) 
is not consistent with the assumptions we made to derive Eq. ( j3.l| ), and therefore it 
does not approximate a physically existing stationary front. If _F > F^, the only 
stationary solution of Eq. (3.1) is the unphysical one. The critical field F^ is such 
that the expansion of the right hand side of ( |3.l| ) about the two coalescing stationary 
solutions has zero linear term, 2 + Ag'{uQ) — 0, and 



(3.2) 



2uo + Ag{uo) - Ui 



Us 



Fc. 



These equations for Fc and uo{A, Fc) have been solved for the FK potential, for which 
uo = cos~^{—2/A) and Ui + U3 = 2sin~^(Fc/A) + 27r. The results are depicted in 
Fig. 3.2, and show excellent agreement with those of direct numerical simulations for 
A > 10. Our approximation performs less well for smaller A, and it breaks down at 
A = 2 with the wrong prediction Fc = 0. Notice that Fc{A)/A ^ 1 as ^ increases. In 
practice, only steady solutions are observed for very large A. 

Let us now construct the profile of the traveling wave fronts after depinning, for 



^^a{F-Fc)+Pvl 



F sligthly above Fc- Then uo(t) 
(3.3) 
(3.4) 

(3.5) /3 



Mo (A, Fc) + vo{t) obeys the following equation: 



a — 1 



1 



Ag'{U,{Fc/A)) Ag'{Us{Fc/A)) ' 



where we have used 2 + Ag'{uo) — 0, (3.2) and ignored terms of order {F — Fc)vo 
and higher. These terms are negligible after rescaling vq = {F — Fc)^ip and r = 
(F ~ Fc)^t. The coefficients a and /? are positive because g'{Ui) > for i = 1,3 and 
g"{uo) < since uo S (f/i(0), J72(0)). For the FK potential, a = l + 2/y/A'^ - F^ and 
P = _ ^2/2. Equation (Q has the (outer) solution 



(3.6) 



vn{t) 



la{F~Fc) 



tan 



(Va/? {F^Fc){t-h)) 



which is very small most of the time, but it blows up when the argument of the 
tangent fimction approaches ±7r/2. Thus the outer approximation holds over a time 
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interval {t - to) 7r/^a/3 {F - Fc), which equals ny^^TjaiA^ - 4)- i {F - Fc)^ 2 for 
the FK potential. The reciprocal of this time interval yields an approximation for the 
wave front velocity, 



(3.7) 



ciAF) 



^aP{F~F,) 



or c ( - 4) 3 ( 1 + 2/ „ ^2) I (i^ - i^J I / (tt V2) for a FK potential. The minus 

sign reminds us that wave fronts move towards the left for F > Fc- In Figs. 3.2(a) 



and (b) we compare this approximation with the numerically computed velocity for 
A = 100 and A = 10. 



back to (3.1 
Fc and wo(i) 



When the solution begins to blow up, the outer solution (3.6) is no longer a good 
approximation, for uo(i) departs from the stationary value Uf){A,Fc). We must go 
and obtain an inner approximation to this equation. As F is close to 
— uq{A,Fc) is of order 1, we solve numerically (3.1) at = i^c with 
the matching condition that uo{t) - uq{A, Fc) ~ 2/['k.J l3/[a{F ~ Fc)] - 2/3 (i - to)], 
as {t — to) — > —00. This inner solution describes the jump of wq from uo{A,Fc) to 
values on the largest stationary solution of Eq. (3.1), which is close to During 
this jump, the motion of uo forces the other points to move. Thus, U-i{t) can be 
calculated by using the inner solution in (1.1) for uo, with F = Fc and m_2 ~ f/i- A 
composite expansion Q constructed with these inner and outer solutions is compared 
to the result of direct numerical simulations in Fig. 3.4. 



Notice that (3.3) is the normal form associated with a saddle-node bifurcation 
in a one dimensional phase space. The wave front depinning transition is a global 
bifurcation with generic features: each individual point Un{t) spends a long time, 
which scales as \F — i^cl"^, near discrete values u„(A, Fc), and then jumps to the 
next discrete value on a time scale of order 1. The traveling wave ceases to exist for 

F < Fr. 
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3.2. Theory with several active points. The approximations to Fc{A) and 
the wave front speed provided by the previous asymptotic theory break down for small 
A. In particular, for the FK potential and A < 2, no double zeroes of 2a:: + ^ sin(a;) — 
{F + Ui + U3) are found for F — Fc- What happens is that we need more than one 
point to approximate wave front motion. Depinning is then described by a reduced 
system of more than one degree of freedom corresponding to active points. There is 
a saddle-node bifurcation in this reduced system whose normal form is of the same 
type as ( ^.3[ ). The jump of the active points after blow up is found by solving the 
reduced system with a matching condition JlOt . 

We explain our procedure for a finite number of active points. The front is 
formed by two tails, very close to Ui{F/A) and {/^{F/A) respectively, and several 
intermediate points . . . ,um- The reduced system describing the dynamics of the 
front for F > Fc is: 

dui 

— = u,+i - 2uj + - Ag(Ui) + F, -L, . . . , M, 

(3.8) u^L-i = Ui{F/A), UM+i - U3{F/A) 

At F = Fc this system has a stationary solution ...um- Writing Ui{t) = Ui + Vi{t), 
we obtain for Vi, 

dv ■ A 

(3.9) = Vi+i - 2vi + v^^i - Ag'{ui)Ui - —g"{ui)vf + F - F^ i = -L, ...,M 

F- Fc F~Fc 



Ag'{UiiFc/A)) ' Ag'iUsiFc/A)) 
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(a) A=100,N=1 
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Fig. 3.4. Comparison of asymptotic and numerically calculated wave front profiles near Fc: 
(a) Complete wave front profile as indicated by the trajectory uo{t). (b) Zoom near the largest jump 
in the profile, (c) Zoom near the jump preceding the largest one after translating the asymptotic 
profile. The latter has been calculated by inserting the approximate uo{t) in the equation for U—\{t) . 



The tridiagonal matrix Ai defined as 

(2 + Ag'{u^L) -1 

-1 2 + Ag'{ui^L) -1 



(3.10) 



V 



-1 2 + Ag'{uM-i) -1 

-1 2 + Ag'{uMy 



has L + M strictly positive eigenvalues plus a smallest eigenvalue A ~ with an 
associate positive eigenfunction V . We choose V such that Y^^L-l ~ ^ ^^'^ write 
an eigenfunction expansion for v: 



M 



(3.11) 



v{t) ^ LpV + ^ Wiipieyi^{-\it). 



Thus v(t) ~ Vip{t) as time increases. Let D = — -j diag [g" {u^l), . . . ,g" {um)), 
F2 = (K_\,...,y|,)and 

w = {F- F,) (l 



Ag'{U^{F,IA)) 



,1,...,1,1 



Ag'{U^{F,/A)) 



Then system ( p.9D becomes 

dip 



(3.12) 



dt 



V - ifMV + Lp^ DV^ +w^if^ DV^ + w, 
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as time increases. Multiplying by the transpose of V, we get an evolution equation 
for the amplitude (p{t): 



(3.13) 

where now a and /3 are: 

M 



i=-L 



V-L , Vm 



Ag'iU.iFjA)) Ag'iU^iFjA)) 



>0, 



A 



M 



i=-L 



The coefScient a is positive because g'{Ui) > for i = 1, 3. We have checked numeri- 
cally that /9 > for different nonlinearities and values of A. An intuitive explanation 
follows. First of all, notice that g"{u) > for u e (C/2(0), ?73(0)) and g"{u) < for 
u e {Ui{0), 1/2(0)) ■ For large A, the largest component is Vq, the others are negligi- 



ble and we have one active point as in the previous subsection; see Fig. 3.4. Then 
P ~ -g"{uoiA,Fc))VQ > because uo < C/2(0), which implies .g"(uo) < 0. As A 
decreases, Vo is still the largest component and g"{uo{A, Fc)) < 0. Now there may be 
other terms with g"{ui{A, Fc)) > and we only have numerical evidence that /3 > 0, 
but not a proof. 



Notice that (3.13) is the normal form of a saddle- node bifurcation. Its solution is 
again (3^), which blows up at times {t — to) — ±l/(2c), where 



(3.14) 



c(A,F) 



— ^aPiF-Fc), 

TT 



as discussed before, c is the wave front speed near Fc, approximately given by the 
reciprocal of the time during which the outer solution holds. 

Figures 2.1, 3.3 and |3.5| show the critical field, wave front velocities and profiles 
for different values of A G (1, 10) corresponding to the FK model. We have compared 
results of direct numerical simulations to those of our theory for = L -|- M -I- 1 active 
points. Provided N — L + M+1 active points have been selected, we find the smallest 
eigenvalue of the matrix Ai and move F until X{F, A; N) = 0, N ^ L + M + 1, which 
yields an approximation for Fc{A). See Fig. ^.1[ . The wave front velocities can be 
calculated by means of Eq. ( p. 14 ) and have been depicted in Fig. 3.3. 

The wave front profiles near Fc can be determined as follows. We start with an 
initial condition, m„(0) « UniA, Fc), or 1^9(0) = in Eq. ( p.ll ). The active points blow 
up a.t t ^ ±(2c)~^, for example as 



(3.15) 



,{t)^Ur,{A,Fc) + 



V, 



provided t ±l/(2c). At these times, we should insert a fast stage during which the 
Un{t) are no longer close to u„{A,Fc), as an inner layer. The inner layer variables 
u„(t) obey Eq. (1.1) with F = Fc and the boundary conditions Un{t) Un{A, Fc) 
[according to Eq. (3.15)] as i — > —00 and u„(i) Un+i{A, Fc) as t ^ 00. To get 
a uniform approximation, we notice that the blow up times are tm — (2c) +171/0, 
m IE Tj. Let us denote by (t), t ~ (t — tm), the solution of Eq. (1.1) with 
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F — Fc and the boundary conditions ul™'' (r) 
,(™) 



^{A,Fc) as T 



-oo and 



ufi (t) ^ ■u„+m+i(yl, Fc) as T — > oo. During the time interval {t-L-n-i,tM-n) = 
(-(2c)-i - (L + n)/c, (2c)-i + {M -n)/c), that u„(i) needs to go from Ui{Fc/A) to 
U3{Fc/A), the uniform approximation to the wave front is 



M-n 



(3.16) un{t) 



m—~n 



^ {u(™)(t-i„)+«(r-i)(t-t„,_i) 

L-l 



(A^^c) 



V 



Then u„(t_L_„_i) - L/i(FcM) and u„{tM-n) - UsiFjA). In Eq. ( |3.16D , the indica- 
tor function X(t„_i,t„,) is 1 if tm-i <t<tm and otherwise. Therefore, X(tm_i,tm) = 
0(tm — i) ~ S{t„i^i — t), where (?(a;) = 1 if a; > and otherwise. Written in terms 
of the variable z ~ n — ct, such that Un{t) — u{z), u{z) = Un{(ri — z)/c) = uq{—z/c). 
Then Eq. (3.16) becomes 



M . 

(3.17) u{z) ^ 

m=-L-l ^ 



z + m + 



c 

Z + 771 



(rn-l) 



Z + 777 - 



/3(z + r7i + i) f3{z + m-^) 

















e 


(- 


- 777 - 


-0 


- e 




h 777 -1) 



for -M - 1/2 < z < L + 1/2. We have u{L + 1/2 + 0) ~ Ui{F^/A) and u{-M - 
1/2 — 0) ^ J73(Fc/^), and therefore Eq. ( 3.17 ) approximates the wave front profile. 
In Fig. |3.5| , we have depicted the wave front profile in two ways, by drawing uo{t) and 
u{z) = uq{—z/c). Notice that the largest source of discrepancy between numerical 
calculations and our asymptotic approximation is the error in determining the wave 
speed. The discrepancies are more evident for Uf){t) because of the different horizontal 
scale used to depict u{z). 

How do we determine the optimal number of active points? For large enough 
N = L + M + 1 and a given A, the eigenvector V corresponding to the smallest 



eigenvalue of the matrix M in (3.10) has a certain number of components that are of 
order one, whereas all others are very small. The number of components of normal 
size determines the optimal number of active points: only one point if A is larger 
than 10, five if A = 2, etc. Keeping less active points than the optimal number 
results in larger errors whereas keeping more active points than optimal does not 
result in a significantly better approximation. The eigenvector of the reduced system 
of equations for the active points is a good approximation to the large components of 
the eigenvector corresponding to the complete system. As we approach the continuum 
limit, more and more points enter the reduced system of equations and exponential 
asymptotic methods become a viable alternative to our methods. 

3.3. Depinning transition as a global bifurcation. We have shown that 
the depinning transition is a global bifurcation in a reduced system of equations 
corresponding to the active points. Starting from a stable stationary solution, the 
smallest eigenvalue of the system linearized about the stationary solution becomes 
zero at the (approximate) critical field, and its associated eigenfunction is positive. 
The stationary solution disappears as the critical field is surpassed. Beyond it, the 
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Fig. 3.5. Comparison of theoretically predicted and numerically calculated wave fronts near 
Fc for A = 2 using N = 8 active points: (a) Trajectory of one point, (b) wave front profile, 
u{z) = uo{^/|c|). 



active points M„(i) spend a long time, of order {F — F(.)~^ , near the stationary values 
u„(A, i^c), and then jump to u„+i(A, i^c) on an order 1 time scale. Near the critical 
field, the depinning transition is described locally by the normal form of a saddle-node 
bifurcation. For F > Fc (or F < —Fc), the bifurcation amplitude blows up in finite 
time, on a time scale of order | — Fc|^2. The construction of the wave front profile 
is completed matching the outer solution given by the saddle-node normal form to a 
solution of the reduced system of active points at = i^c (a mathematically related 
phenomenon occurs in a mean-field model of sliding charge-density waves Q ) . 

We conjecture that the depinning transition in the infinite system (lA) is a global 
bifurcation of the same type as for the reduced system of active points. At the critical 
field, two stationary solutions of Eq.(^]^) (one stable, the other unstable) coalesce 
and disappear. For F > Fc (or F < —Fc), the wave front profile is constructed as 
indicated above for the reduced system. To prove this conjecture, we could repeat our 
construction in Section 3.2 for an infinite number of points. This is possible because we 
know that the infinite system, linearized about the 'stable' steady solution Ui{A, Fc) 
at F — Fc, has a zero eigenvalue and an associated positive exponentially decaying 



eigenfunction V. Using V, we obtain the normal form equation (3.13), where now 



(3.18) 



a 



^ l^, > 0, (3 



A 
" 2" 



OO 

E 

i— — oo 



We should now prove that the coefficient (3 is positive and that the infinite system 
has solutions connecting Un{A, Fc) to Un+i{A, Fc) and satisfying the matching con- 
dition. We justified that /3 > for the finite system in Subsection ||.2 and we show 
in Proposition A.l (Appendix A) that the eigenfunction for the infinite system can 
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be approximated by the corresponding eigenfunction of the reduced system with a 
finite number of active points. The existence of travehng wave solutions for F > 
ensures that the infinite system has solutions connecting Un{A, Fc) to Un+i{A,Fc). 



The velocity of a wave front in the infinite system is again given by (3.14) with the 



coefficients (|3.18| ) 



Remark 5. By using comparison techniques, it is possible to prove that solutions of 
discrete RD equations with finitely many points and Dirichlet boundary conditions 
approximate solutions of the same equations with infinitely many points. In the 
continuum limit, the wave fronts approach constant values exponentially fast as i — > 
±cxD. This exponential decay justifies the active point approximation in two ways. 
First of all, the number of active points needed to approximate well the wave fronts 
of the infinite system decreases as A increases. It is usually better to add another 
active point to the approximate system than to patch rigid tails to the last active 
points of a wave front by generalizing Kladko et al's active site approximation [ p8| . 
Secondly, exponential decay at the ends of a wave front causes the operator of the 
linearized problem about the wave front to be compact, and therefore to have discrete 
spectrum (see Appendix A). This fact justifies that the normal form we calculate by 
using active points approximates the correct local normal form of the depinning global 
bifurcation. 

4. Conclusions. In this paper, we have studied depinning of wave fronts in dis- 
crete RD equations. Pinned (stationary) and traveling wave fronts cannot coexist for 
the same value of the forcing term. There are two different depinning transitions, i.e. 
two different ways in which a pinned front may start moving. The normal depinning 
transition can be viewed as a loss of continuity of traveling front profiles as the critical 
field is approached: below the critical field, the fronts become pinned stationary pro- 
files with discontinuous jumps at discrete values u„. The wave front velocity scales as 
\F — Fc| 2 near the critical field Fc- For sufficiently large A (far from the continuum 
limit), the critical field and these fronts can be approximated by singular perturbation 
methods which show excellent agreement with numerical simulations. These methods 
are based upon the fact that the wave front motion can be described by a reduced 
system of equations corresponding to the dynamics of a finite number of points only: 
the active points. 

Besides the normal depinning transition, certain nonlinearities present anomalous 
pinning (pinning failure): the velocity of the wave fronts is not zero except at zero 
forcing, just as for continuous RD equations. These nonlinearities are characterized 
by smooth profiles of stationary and moving wave fronts, by having zero critical field 
and by a linear scaling of wave front velocity with field. 

Appendix A. Characterization of the depinning threshold. 

In this section we establish the 'depinning criterion' which provides a characteri- 
zation of Fc{A): 

Theorem A.l Set F ^ and A > 0. Assume that the nonlinearity g £ has 
three zeroes Ui, Ui < U2 < U3, is odd abou t U2 and satisfies g'{Ui) — g' {U^) > 0. 
Let Un he a stationary increasing solution of ^l.ij ), symmetric about U2 and such that 
it_oo — Ui and Uqo — U3. Let Ai(A,0) he the smallest eigenvalue of the zero field 
operator Cq of at F ^0: 

- {Vn+l - 2u„ + Vn-l) + Ag'{Un)v,i = Xi{A, 0)Vn 

(A.l) v±n exponentially as n — > cx) 
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If Xi{A,0) > 0, then FJA) > and for \F\ < Fi.{A) there exist increasing stationary 
solutions Un{A, F) of ([/jj with u^rx- — Ui{F/A) and — U^{F/A). Moreover, the 



smallest eigenvalues of the operator C{F), corresponding to the linearization of (1.1) 
about Un{A, F), 

- {Vn+l - 2Vn + Vn-l) + Ag' (UniA, F))Vn = \l{A, F)Vn 

(A. 2) v±n — > exponentially as n ^ oo, 

are strictly positive for \F\ < Fc{A). We can characterize Fc{A) as the zero of the 
smallest eigenvalue, \i{A,Fc{A)) — 0. 

This theorem wiU be proven in subsection A. 2. T o cal culate \i{A, F) and Fc{A), 
we approximate the infinite tridiagonal matrix in Eq. ([A.2| ) by a x A'^ matrix, where 
A^ is the number of active points. Similar truncation approximations were used in ||2^ 
to calculate the lowest eigenvalue of an infinite tridiagonal matrix. For values of A 
which are not too small, numerical simulations show that the matrices M in (|3T0| ) 



have positive eigenvalues. The eigenvector V{A,F,N), (chosen to have norm 1) cor- 
responding to the smallest eigenvalue, X{A,F,N), is positive and it is 'concentrated' 
in the central components, V^_m(^), . . . , V^(yi)- AH other components are very small. 
The number of significant components m{A) does not change as N increases, but it 
increases as A decreases. For large A, m{A) — and only Vq is significant. Pro- 
vided A^ is large enough, the eigenvalues \{A, F, A) and the eigenvectors V{A, F, N) 
approximate well the smallest eigenvalue and associated eigenfunction of the infinite 
problem, as we indicate in the next subsection. For a fixed value of A^, the eigenvalues 
X{A, F, N) decrease as A decreases. For fixed A^ and A, they decrease as F increases 
from F = to values close to Fc{A). In the next subsection, we collect several results 
on eigenvalues for this type of problems. 

A.l. Eigenvalue problems. Before proving Theorem A.l, we should make sure 
that our linear operators do have eigenfunctions and eigenvalues. We consider the real 
valued and symmetric operators C{F) u„ = Ag'{un)vn — {vn+i — 2wn + w„_i), in spaces 
of sequences decaying exponentially at infinity. Their spectra are discrete and real 
(these operators are compact), and we would like to make sure that they are not empty. 
Since we are interested mainly in the smallest eigenvalue, we shall use its variational 
characterization, prove that this eigenvalue exists and characterize its dependence on 
the parameters A and F. We shall also describe finite dimensional approximations of 
eigenvalues and eigenfunctions. 

Let us first look for necessary conditions for X{A, F) to exist. Let A G II be 
an eigenvalue of C{F) with eigenfunction Vn- Multiplying L{F) Vn — Al/„ by and 
summing over n, we obtain 

= 5](K+i - Vr.f + {Ag\u^) ~ \\Vl 

n 

(A.3) > ^(K+i - Vnf + [Amin„(g'(ti„)] - A) ^ V^. 

71 71 

Thus, A = - K)' + Ag'{un)V^]/Y.V^ > Amin„g'(u„). This inequality 

implies that A is positive if u„ does not take on values in the region where g' is 
negative, which occurs for large enough A > 0. In general, we can only say that 
A > Ag'{U2) for g' attains its minimum value in [Ui, C/3] at U2- Similarly, 



= ^(K+i - Vn f + [Ag'iun) - \]V^ 
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(A.4) < [4 + Amax„(5'(M„)) - A] ^ 

n 

Therefore, A < Ag'{Ui) + 4 = Ag'{U3) + 4, for g'{u) attains its maximum value in 
[Ui, U3] at the end points, u ~ Ui and u — U3. 

The smallest eigenvalue A is given by the Rayleigh formula 

(A.5) A ^ Min EJ-'n+.-w^r + Ag'Mwl 

En K 

where the infimum is taken over our space of exponentially decaying functions. That 
Minimum is attained at an eigenfunction Vn solving (A.l). Now, (A.l) may have 
solutions decaying at ±00 only if the difference equation 

(A.6) Vn+i + (-2 - Ag'iUi) + A) K + K-i = 

has solutions of the form r" with r < 1. This happens when (—2 — Ag'{Ui) + A)^ > 4. 
Thus either A > 4 + Ag'{Ui) > (excluded above) or A < Ag'{Ui). We conclude 
that A < Ag'{Ui) is a necessary condition to attain the Minimum (A. 5) at a positive 
eigenfunction decaying exponentially at infinity. 

We now establish sufficient conditions for the Minimum (|A.5| ) to exist. 

Lemma A.l.- Conditions for the existence of positive decaying eigen- 
functions Let F ~ 0, A > 0, and let the nonlinearity g satisfy the hypotheses in The- 
orem A.l Let Un be a stationary increasing solution of jl-j ) such that 7i_oo — Ui{0) 
and Uao — £^3(0). Given an exponentially decaying sequence w = Wn, we define 

J(u;) ^ T.nii^n +1 - Wnf + Ag' {Un)wl] 



Let us suppose that there is a sequence Wn such that J{wn) < Ag'{Ui). Then the 
infimum 

f.rj^ \ Tf Y.ni'"n+1 - Vnf + Ag' {U^)vl] 

(A-7) A = Inf^^j^-.„,,.<^ 

is attained at a positive function Vn which decays as r(yl, A)'"' at infinity, with < 
r(AA) - [2 + Ag'{Ui) - A - ^{-2 - Ag'{U^) + A)^ - 4]/2 < tq < 1. A = Ai(AO) 
and Vn solves ( |^.4 J. 

Remark A.L The value < ro < 1 is determined in the proof. Note that r{A, A) is a 
decreasing function of A but an increasing function of A. 



Remark A. 2. We have shown above that Ag'{U2) < A < Ag'{Ui). Thus the smallest 
eigenvalue shrinks to zero as A ^ 0, although we do not have a proof that it does so 
monotonically. 

Proof: Clearly J{w) is bounded from below by Amin„((7'(u„)). We choose 

ro = r{A,J{wn)) € (0,1) and define ||w||o = E^o'^'"'k"P- Let w"' = w^f , m > 
be a sequence minimizing J{w): \\w'^\\q < 00 and J{w"^) — > A when m 00. We 
replace with ^ v^^ = w^/Ww'^Wo. Then, \\v"^\\a = 1 and J(u") J(w") ^ A. 
||v™||o = 1 imphes that, uniformly in m, < rj,"' and J2n>n{e) l^rTP < ^ fo'' 

n(e) large enough. Thus, a subsequence tends to some limit V — Vn such that 
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\Vn\ < Tq"' and J2n l'^™ ~ tends to zero as m tends to infinity. Therefore, 
J(u™) — > J{V) = A and the infimum is attained at the sequence V = Vn- Moreover, 



V = Vn satifies the Euler equation (A.l) for the minimization problem, which then 
imphes that Vn decays as stated in the Lemma. 

On the other hand, J(|Ki|) < J{Vn) and we can choose nonnegative Vn- But then 
A has to be the smallest eigenvalue Xi{A, 0). 
■ 

Lemma A. 2.- Choice of the sequence w„ with J{wn) < Ag'{Ui). Let F — 0, 

A > 0, and Un he as in Lemma A.l. Let u{x) be the solution of the boundary value 
problem d^u/dx^ — g{u), with u{-~oo) — Ui, u{oo) = U3 such that u{0) = U^- 
• For sufficiently small A < \, we have 

(A. 8) \un+i - Un\ max[„^^„^^^] \g"\ < 2g'{Ui), Vn, 

and Wn — Un+i — Un satisfies J{wn) < Ag'{Ui). An estimation of the appro- 
priate values of A indicates that they should be smaller than 



A< 



max[u^^u^]\g"\max^\du/dx\ 



• For A > 2g"'{Ui), we can choose Wn = for \n\ > M > 0, Wn = for 
n > and Wn — r"+^ for n < 0. 
Proof: The function Wn = Un+i — Un is a solution of 

O , .3(Wn+l) - 5(Mn) 

W„+x - 2Wn + Wn-1 = A Wn 

Un+1 ~ Un 

that decays at infinity as r{A, O)'"'. Multiplying this equation by Wn and adding over 
n, we obtain 

(A.9) > {[Wn+l-Wn) +A w„]^0. 

^ \ - Un J 

This result can be used to calculate J(w„): 



(A.IO) JiWn)=A- 



by the mean value theorem. Thus, J{wn) < {A/2) max„ (\un+i—Un\ inax[u„,u„+i]\g" \) ■ 
For sufficiently small A, |u„+i — u„| < C^/A, so that J(wn) < CAi uYax\g"\/2 < 
Ag'{Ui). More precisely, for small A, w„ — w„+i — w„ ~ u{{n + 1)^/A) — u{n^/A) ~ 
v^w'(^). Then |m„| < max \du/dx\y/A. 

To prove the other case, we observe that J(w„) = (r — 1)^ + ^Ei^o 5''("n)''^" 
smaller than ^^'(ui) provided (l-r)/(l + r) < Ag"'{Ui)/[2{l - {r{A,0)r)^)] with 
r{A, A) defined as in Lemma A.l. The last inequality holds ii A > 2g"'{Ui). 



Remark A. 3. For the FK nonlinearity, the first condition of the Lemma holds for 
A < 0.9, and the second condition for A > 2. For intermediate values, numerical 
simulations show that w„ — Un+i — u„ satisfies J{wn) < Ag'{Ui). 
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Proposition A.l . - Finite dimensional approximations Let Un{A,F) be a 
stationary solution of jj.j j under the hypotheses in Theorem A.l for \F\ < FdA). Let 
Ai(A, i^) the smallest eigenvalue of the operator C{F) (linearized about Un{A, F) ) and 
X{A, F, N) be the smallest eigenvalues of the matrices ( 3.1(\ ). Then, X{A,F,N) 



Xi{A, F) as N oo. As a consequence, if V > is an eigenfunction associated to 
Xi{A,F) with J2n^n — 1 '^'^'^ ^i^) > are eigenvectors associated to X{A,F,N) 
such that J^n VniN)^ = 1, V{N) -^V as N ^ oo. 

Proof: It follows from the Rayleigh characterizations for the smallest eigenvalues: 

E!°oo [K+i - ^n? + Ag'{u^{A, F))wl] 



(A.ll) Ai(A,F)=Mmy^_,|„|^, =^3o 5 

[K+i - Vn? + Ag'iuJA, F))v^] 
(A.12) XiA, F, N) = Min ^^^^^ ^ 'J y \ ny , >> n\ 

Letting Wn ~ Vn for n = —L, . . . , M and w„ = otherwise, we see that Xi{A, F) < 
X{A, F, TV), N ^ L + M + 1. Let now w„ be an eigenfunction for Ai(A, F) such that 
E^^wl^l. Then, 

^ m ^ ~ E„<-L,„>A./[K+i - w,,)^+Ag'{u^{A,F))wl] 
X[A,l',l\) < ^:^jg . 



L 



We conclude that X{A, F, N) Ai {A, F) as ^ oo. This and the exponential decay 

of V prove the convergence of the eigenvectors. 

■ 

A. 2. Proof of Theorem A.l. The theorem will be proved in two steps and 
for simplicity, in the particular case of periodic g. In this case, U3{F/A) — C/3(0) = 
Ui{F/A) — Ui{0), which allows us to use symmetric sub and supersolutions. Small 
modifications are required in the general case. 

Step 1: Fc{A) > 0. 

We use the existence of a positive eigenfunction w„ associated to a positive eigen- 



value Xi{A, 0) to construct stationary supersolutions for (LI) when F > small. The 
known solution u„ provides a stationary subsolution. 
We look for a supersolution of the form 

(A.13) vjn^un + {l + S){Ui{F/A) - f/i(0)) + ew„ 

with (5 > to be chosen and e, F small to be determined. Let us check that 

(A.14) Wn+l ~ 2Wn+Wn^l < g{Wn) - F, Woo > J/g (F/A) , W_oo > Ui{F/A) 

holds. The conditions at infinity are satisfied for any S > 0. Provided 
(A.15) 1(1 + S)(Ui{F/A) - C/i(0))| < ke, 



inequality ( A.14 ) holds if 
e [vn+i - 2vn + Vn-i) < Ag' (un) [ev„ + (1 + 6){UiiF/A) - C/i(0))] -F + OiAe^). 



Using (A.l) we are left with 

F < eAi(A,0)f„ + ^(1 + 6)g'{un){UiiF/A) - C/i(0)) 
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Now, Ui{F/A) = g ^(^)' ^^'^ inverse being taken near C/i(0), in the region with 
g' > 0. Using g^^{x) ^ g^^{xa) + {g"^)' {xq){x — Xq) we obtain 

(A.16) UiiF/A)^g-\F/A)^UiiA,0) ^ 



Ag'iUiiO)) ■ 

Thus, the condition for w„ to be a supersolution is 

F < eXi{A,0)vn + g'iun)-^j^^^^^F 

Let M be sufficiently large. We distinguish two different ranges of indices n: 

• For \n\ > M, g'{un) > and w„ <C 1. Then the right hand side of the previous 
inequality is dominated by the second term. We choose 5 large enough to 
ensure F < g' {un){l + 5) F / g' {Ui{Q)), that is, f + ,5 > g' {Ui{0))/ g' {um). 

• For small g'{un) < 0. The previous inequality is satisfied provided we 
choose F so small that 

1 + I.9'K0I^7^7^) ^<eAi(A,OK, 
for a fixed value of 5. 

With these choices, w„ satisfies ( A.14|) . Note that these choices are compatible with 
condition ( |A.f5| ). Using ( |A.16| ), ( |A.f5| ) becomes [1 + 5) F /[Ag' {Ui{0))\ < ke. This 
holds for small enough F. 

Let F > be small enough for a w„ defined in (A. 13) to be a supersolution with 
S,e adequately selected. Now, let hn{t) be a solution to (1.1) for such F > with 
initial datum hn{0) satisfying u„ < hn{0) < Wn- Then, u„ < /i„(i) < w„, for all 
i > 0. Therefore, propagation is excluded and the solutions are pinned. 

Stationary solutions Un{A, F) for such F > can be obtained as long time limits 
of solutions hn{t) to ( pTl] ) when /i„(0) is increasing, tends exponentially to Ui{A, F) 
(resp. U3{A,F)) at — oo (resp. oo) and u„ < ft-„(0) < w„. We conclude that Fc{A) > 
0. 

Step 2: Xi{A,F) > for |F| < Fc{A) and Xi{A, FdA)) = 

To fix ideas, we take F > 0. The case F < follows by symmetry. From Step 1, 

we know that Fc{A) > and there are stationary solutions Un{A, F) of (IT) exist for 

F > small, that are increasing from Ui{F/A) to U3{F/A). 
In an analogous way as we did for F = we get: 



(A. 17) Xi{A,F) = Min 



^„[(w„+i - u;„)2 + Ag'{un{A, F))wl 



^0 w^<Cqo 



w, 



2 



This formula defines Ai [A, F) as a continuous function of F. That Ai (A, 0) > implies 
Ai(A, F) > up to some F^, at which Ai(A, F^) = 0. As long as Ai(A, Fi) > 0, we can 
obtain stationary solutions for F > Fi (close to Fi ) , as done in Step 1 . This procedure 
cannot continue forever since such stationary solutions do not longer exist for F close 
to A: eventually g{U) — F/A ceases to have three solutions and the stationary wave 
fronts cannot be constructed. Thus, we must reach a value Fc at which Xi{A, Fc) = 0. 
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